Wave Localization and Density Bunching in Pair Ion Plasmas 
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By investigating the nonlinear propagation of high intensity electromagnetic (EM) waves in a pair 
ion plasma, whose symmetry is broken via contamination by a small fraction of high mass immobile 
ions, it is shown that this new and interesting state of (laboratory created) matter is capable of 
supporting structures that strongly localize and bunch the EM radiation with density excess in 
the region of localization. Testing of this prediction in controlled laboratory experiments can lend 
credence, inter aha, to conjectures on structure formation (via the same mechanism) in the MEV 
era of the early universe. 
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The physics of pair plasmas was turned into an even 
more exciting field of investigation when it descended 
from its astrophysical heights to the terrestrial labo- 
ratory. The story of laboratory pair plasmas had its 
defining moment with the successful creation of a " suffi- 
ciently" dense pair-ion (pi) plasma - consisting of equal- 
mass, positive and negative fuUerene ions (C^ and Cqq) 
Unlike the electron-positron (e-p) plasma systems 
(both of the astrophysical 0, H, 0, [1, 0, 13 laboratory 
[Sj variety), the fuUerene plasma has a long enough life 
time that the collective behavior peculiar to the plasma 
state can be experimentally investigated under controlled 
conditions. 

The dynamics of pair plasmas is expected to be dif- 
ferent from the standard electron-ion plasma where the 
different mass of the species, automatically breaks the 
symmetry between the constituents. On the other hand, 
if the pair plasmas are prepared in identical conditions, 
they must remain symmetric - for example their thermal 
speeds and temperatures are likely to be similar. 

This is, indeed, the case as evident from the tremen- 
dous surge in theoretical activity [1, [13, EH 
to interpret and understand the experimental results [ij. 
Most of these papers attempt to explain the experimental 
findings in terms of the linear and nonlinear properties 
that may be accessible to only pure pi plasmas. Some- 
what different problem appears to be when in PI plasma 
while its creation a significant fraction of free electrons 
remain and specific linear modes may develop in such 
plasmas [iTI. 

The aim of this paper, however, is totally different. 
We wish to propose, here, experiments that explore phe- 
nomenon that are associated with slightly contaminated 
plasmas - where the symmetry is broken, say by the pres- 
ence of a small amount of a heavier ions (heavier than 
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the ions that constitute the main plasma). Our motiva- 
tion is two fold: 1) to study, per say, a possible nonlinear 
bunching of electromagnetic waves through this new and 
exciting state of matter whose composition may be highly 
controllable (naturally at some future date), and 2) to 
create (from this restricted perspective) an approximate 
replica of the cosmic plasma of MEV era which consists 
primarily of electrons and positrons with a small concen- 
tration of symmetry-breaking ions ^5, 0) 01 ■ T'^^' hope is 
that the laboratory experiment will show that, in a min- 
imally contaminated pair plasma, electromagnetism can 
provide mechanism for density bunching lending further 
credence to the idea of the electromagnetic origin of the 
large scale structure of the universe. 

It is important to stress that, although cosmic consid- 
erations do provide a motivation, the principal intent of 
this paper is to investigate nonlinear phenomena in this 
potentially very interesting, versatile, and controllable 
new state of matter. Because of relatively low densities 
and heavier mass of the species, the corresponding fre- 
quencies and scale lengths will be in a novel range quite 
different from the c-p plasma counterpart. One of the 
challenges, therefore, will be to derive laboratory condi- 
tions suitable for the observance of bunching of the EM 
waves. Similar behavior could be expected in doped (or 
dust-contaminated) fuUerene plasmas in laboratory" [7|. 

The fuUerene plasmas of massive ions, however, are not 
quite suitable for this experiment since all the frequencies 
associated with the collective modes (plasma frequency, 
acoustic and Alfven frequencies) tend to be rather low. 
Fortunately, the group of Hatekayama and Oohara have 
already made considerable progress in the production of 
the hydrogen, H'^-H~ plasmas [H, [l^. Since the initial 
report, both the quality and quantity of this light pair ion 
plasma has been steadily improving. We will not dwell 
on how the H^-H~ plasmas are or will be produced and 
diagnosed (see [i^[ia|). We will simply assume that, in 
a not too distant future, it will become possible to "dial" 
in an H^-H~ plasma with densities sufficiently "high" 
that collective modes can be experimentally excited and 



2 



studied [17} . 

The main thrust of this paper is to investigate the non- 
Unear interactions of EM waves with a primarily H^-H^ 
plasma containing a small impurity of high mass ions 
(positive or negative). Wc will show that the symme- 
try breaking induced by the immobile heavy ions (the 
main components, and -ff^, will have slightly differ- 
ent ambient densities to insure charge neutrality) leads 
to a finite electrostatic field which, in turn, makes the 
plasma capable of supporting stable localized EM wave 
structures of finite amplitude. 

The pair ion plasma with a heavy ion contamination 
will be non-relativistic both in velocity and temperature 
as opposed to the cosmic plasma it is supposed to mimic 
Although our eventual interest will be in the {H^- 
H~ + heavy ion) plasma, we will analyze an arbitrary 
pair plasma whose principal constituents have a mass m, 
and charge ±q. The equations of motion for the dynamic 
species are: 



E+ i(u± X B) 

c 



(1) 



where q = Z \e\ ; — mu^ (u^) are the hydrodynamic 
momenta (velocities) of the oppositely charged species, 
E and B are the electric and magnetic fields and /dt = 
d/dt + • V are the co-moving derivatives. In Eq. 
([1]), we have invoked a generic equation of state — 
P^(n^) to write the pressure term (l/n^)VP^ as VP^. 

It has been shown earher 0, [1] that: (a) the perfectly 
symmetric system ^ cannot sustain a localized elec- 
trostatic field, and (b) the electrostatic field is essential 
for the localization of an electromagnetic wave passing 
through a pair plasma. Since the creation of localized 
nonlinear electromagnetic structures is the theme of this 
paper, a mechanism for symmetry breaking must be pro- 
vided. We do this by doping the pi plasma with a heavy 
ion impurity (for the electrostatic mode excitations in 
e-p-i plasmas see Q) so that the demands of charge neu- 
trality 



q = qn 



q No. 



(2) 



will cause a difference between the densities of 
the main constituents. Here q' is the charge of the 
heavy, non-dynamic, uniformly distributed, contaminat- 
ing specie, and oo is taken to be the fiducial point where 
the perturbations vanish. Equation ([2]) implies the rela- 
tion qW^ = qn^ + q N^o for ambient densities; it also 
leads to 



n — n = n' 



(3) 



By design, the density iVoo of the contaminant is much 
smaller than the density of the main pair (A^oo ^ 



which, evolve via the continuity equation: 



V • (n±u±) = 



(4) 



The equations of motion, now, must be coupled with 
the Maxwell equations. In terms of the vector (A) and 
the electrostatic {(j)) potentials, the latter may be written 
as (Coulomb gauge V • A = 0): 



c^AA 



c^(V^) = 47rcJ, (5) 



Aip = -47rp, (6) 
where the charge and current densities are defined by: 
p ~ q{n^ — n^) — q Noo', J = g(n+u+ — ?i^u^). (7) 
In terms of the dimcnsionless variables 



u = — , n 

c 



A 



qA 



-r, t — uj^t. 



(8) 

where A = [A; ip; {q~'^T'^)] and uj- = {ATrq"^ n;^/my/'^ 
is the Langmuir frequency of major negative species, the 
defining equations read (the continuity equations retain 
their form ([4])): 

d^A ^ , d , 
^-AA+-(V^)-H 



, 



{l + €)n+ + e] 



~dt 



(9) 



(10) 



(11) 



where the equations of motion ([T]), after standard ma- 
nipulation, have been rewritten in a revealing form that 
contains the generalized flows 



n± = u± ± A; 

the generalized vorticities 

fl± = V X n±; 
and the effective energies 



(12) 
(13) 

(14) 



The symmetry-breaking small parameter e = 
\q Woo I qn'^ <^ 1, appearing in equations ((9))- p0)) . 
will eventually be the source as well as the measure of 
the electrostatic field. 

We are, now, all set to analyze the one-dimensional 
propagation {d/dz ^ 0,d/dx — 0,d/dy = 0) of a circu- 
larly polarized EM wave with a mean frequency Uo and 
a mean wave number ko, 

A± = -(x -I- iy) A{z, t) exp {ikoZ — iojot) -\- c.c. (15) 
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Here x and y are the standard unit vectors, and the 
envelope A{z,t) is a slowly varying function of z and 
t. Using the gauge condition {A^ = in the present 
context), we can calculate 



X rj"*^ = X (z X 



dm 



z u 



dz 



dz 



(16) 



and substitute it in the transverse component of pip to 
derive: 



d 



d 



(17) 



The simplest solution Yl^ = is the most relevant here 
since, for the localized solutions we are seeking, all fields 
must go to zero at infinity. The consequential relation 



(18) 



relates the transverse components of the hydrodynamic 
velocities and the vector potential. 

The longitudinal dynamics is obtained from the z com- 
ponent of pTjl using pB]) . Noting that all terms in this 
dynamics vary on a slow time scale, we can introduce the 
following variables for convenience: ^ = z — Vgt ,t = t , 
where Vg = kg /luq is the group velocity of the EM wave 
packet. Assuming Vgd/d£, ^ d/dr, and integrating, we 
find: 



± ^ ±2 

ut = - ut 



- A\ ± 

2 ^ 



(P±-P±), (19) 



The result holds in the transparent plasma limit {ujq 3> 1) 



for which 



1. In the laboratory experiments, this 



condition, which in dimensional terms demands tJo 3> 
\/2w- {— ujp the plasma frequency) can be easily ar- 
ranged for pi plasmas [l|, [l^, [Ty| . In the non-relativistic 
case under consideration (the velocities are normalized 
to c), the first term on the r.h.s. of is readily ne- 
glected. However the second term j4^= must be 

kept because is allowed to be ^ uf^- In (fTO|) . the 
constants of integration are determined from the bound- 
ary condition that, at infinity, EM fields and the plasma 
momenta vanish. 

Similarly, the continuity equation yields: 



(20) 



From (fTO|) and and the quasi- neutrality condition 

(obtained by neglecting A0 in the Poisson equation (flOl) 
one derives, after straightforward algebra, the following 
major relationships: 



(21) 



and 



2e- 



2(1 -£(/.). (22) 



As expected, the electrostatic field (that will bring in the 
nonlinearity, necessary for localization, in the Maxwell 
equation 0) owes its existence to symmetry breaking, 
and is proportional to the parameter e <C 1. It is clear 
that the electrostatic potential will be absent in pure pair 
plasmas with no contamination. This phenomenon was, 
first, investigated for a relativistic plasmas in a cosmic 
setting 0] . In a pure pair plasma (equal density and equal 
temperature for the two species) the radiation pressure 
imparts equal longitudinal momenta to both the negative 
and positive ions (since their effective masses are equal) 
and, thus, causes no charge separation (n_ = n+ and 
= 0). 

To proceed further, let us determine the transverse 
current defined in ([7]) or ([9|) by resorting to (fT8|) . ((2T|l . 
and (|22p . To the lowest order in the nonlinearity, we find. 



(2- 



(23) 



Notice that the nonlinear term coincides with the non 
relativistic limit of the expression derived in Q . We see, 
that without an asymmetry (e = 0), the nonlinear term 
disappears since the electrostatic potential (responsible 
for this term) </> is now zero. 

When the transverse current, evaluated in ([23|) . is sub- 
stituted into ([5]), we obtain the nonlinear equation for 
the evolution of the slowly varying envelope function A. 
Following the standard methodology (see for exam- 
ple), we can derive the final equation in the parabolic 
approximation: 



2iwo-5- 

OT 



OJ, 



+ -e^ \A\^ A = , 



(24) 



where the wave frequency loq satisfies the dispersion rela- 
tion ujq — (2 -I- e) -f /cq implying — 1 for a transparent 
plasma for which ujq ^ 1 has been assumed. With the 
additional obvious normalization, Eq. (|24p is cast into 
the standard form of the Nonlinear Schrodinger equation 
(NLSE), 



dA d^A 



dr dz 



■e^\A\^A = 0. 



(25) 



It is well known that the NLSE can be solved for soli- 
tons by the inverse scattering method. The usual station- 
ary soliton solution is obtained by letting {Q is a constant 
that corresponds to a nonlinear frequency shift) 



A^AiOe 



(26) 



that leads to the nonlinear ordinary differential equa- 
tion(ODE) 



1 



— 42 



A^O. 



(27) 



4 



With boundary conditions appropriate to a localized so- 
lution, ^ = = dA/dS^ as ^ ^ ±00, Eq. (gT]) yields 

A = A„sech(^^c), (28) 

a lump of the electromagnetic field potential with am- 
plitude Am — \/2/e, and localization width d=^ 
^/2/eAm = 1- In standard measures, therefore, the lo- 
calization width is the collision less skin depth of the 
main plasma {d ^ \skin)- 

We have, thus demonstrated that the pair ion plasmas 
for which the symmetry is broken by a slight contami- 
nation (doping) of a heavier immobile ion can support 
stable localized EM wave structures even in the non rel- 
ativistic limit appropriate to the current and near future 
laboratory experiments. To appreciate the importance of 
theoretical results in the context of a laboratory setting, 
let us invoke a hypothetical H^-H~ plasma of density 
10^*^ cm~^ that is expected to become available in near 
future [13] ■ For these densities, the characteristic 
plasma frequency LOp ~ 1.86 * 10^ corresponding to 
a. Vp = ujp/2TT ~ 30 MHz. The first condition for our 
derivation that the plasma be transparent (luq ^ Wp) can 
be easily satisfied by choosing the EM wave frequency to 
be in the few hundred MHz range. 

Since the localization distance was shown to be about a 
skin depth, for ujp = 1.86* 10* s^^, it comes out to be d ~ 
As/cm = c/up ~ 1.6 m. Thus to experimentally observe 
density and electromagnetic field bunching in a doped 
H^-H~ plasma, the plasma must extend over several 
meters (along the propagation direction). Naturally, for 
higher densities, the localization distance will decrease 
and the length requirement on the experiment will be 
correspondingly relaxed. 

By investigating the nonlinear propagation of high in- 
tensity EM waves in a pair ion plasma contaminated with 



a small fraction of a high mass immobile ions (for sym- 
metry breaking), we have highlighted a very remarkable 
property of this new and interesting state of (laboratory 
created) matter - it can strongly localize the EM radi- 
ation with finite density excess in the region of localiza- 
tion. It must be stressed that this is just one example of 
a vast variety of physical phenomena (many with astro- 
physical consequences) associated with pair plasmas that 
can be explored in controlled experiments. 

It is worth repeating that the equation (|25p obeyed by 
the field envelope A comes out to be the NLSE which 
is known to have exact soliton solutions. Solitons carry 
a large density inhomogeneity, and form non diffr active, 
non dispersive localized structures. Such objects have 
been called "heavy bullets of light", and their possible 
important role in cosmology as a source of structure for- 
mation has been already alluded to 0) 3 ■ The solitons 
are also important for energy transfer investigations in 
plasmas. 

There are several ways in which this simple idealized 
calculation can be extended and augmented. For example 
we could, inter alia: 1) investigate other modes of sym- 
metry breaking, 2) introduce a guide magnetic field (see 
for electrostatic mode envelope excitations in pair-ion 
plasmas doped by massive ions (or dust particles), and 
3) increase the dimensionality by introducing weak trans- 
verse dependence. 
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